Bose-Einstein solitons in highly asymmetric traps 
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We obtain analytic solutions to the Gross-Pitaevskii equa- 
tion with negative scattering length in highly asymmetric 
traps. We find that in these traps the Bose-Einstein con- 
densates behave like quasiparticles and do not expand when 
the trapping in one direction is eliminated. The results can 
be applicable to the control of the motion of Bose-Einstein 
condensates. 



I. INTRODUCTION 

The recent experimental realization of Bose-Einstein 
condensation (BEC) in ultracold atomic gases, QJ^] has 
triggered the theoretical exploration of the properties of 
Bose gases. Specifically there has been a great interest 
in the developement of applications which make use of 
the properties of this new state of matter. Perhaps, the 
recent development of the so-called atom laser || is the 
best example of the interest of these applications. 

The current model used to describe a system with 
a fixed mean number TV of weakly interacting bosons, 
trapped in a parabolic potential V(r) is the following 
Nonlinear Schrodinger Equation (NLSE) (which in this 
context is called the Gross-Pitaevskii equation (GPE)) 
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which is valid when the particle density and temperature 
of the condensate are small enough. Here Uq — Airh 2 a/m 
characterizes the interaction and is defined in terms of the 
ground state scattering length a. The normalization for 
-0 is N — J \tp\ 2 d 3 f, and the trapping potential is given 
by 
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■i n , (rj — x,y,z) being, as usual, constants describ- 
ing the anisotropies of the trap Q . In real experimental 
systems the geometry of the trap imposes the condition 



X x = X y = 1. X z = v z jv is the quotient between the fre- 
quency along the z-direction v z and the radial one v r = v. 
Eq. (|]) is strictly valid in the T = and low density 
limit, but has been validated by different ways for the 
current experimental systems, e.g. by the comparison of 
the experimental j|] an d theoretical low energy excita- 
tion spectra of the condensates |6[0]. Recent theoretical 
work extends the applicability of the GPE to the high 
density limit ||,||. 

When a > the interaction between the particles in 
the condensate is repulsive, as in most current BEC ex- 
periments Jlp|, [lO| , pT| . In opposite case (a < 0), the in- 
teraction is attractive H,|l2|. 

Although the GPE is widely accepted as a valid model 
for the dynamics of the BEC at T ~ K, the knowl- 
edge of the dynamics of the condensates is scarce since 
the GPE is non-integrable and explicit solutions are not 
known. In the positive scattering length case Eq. (jl]) 
has been solved numerically for cylindrically symmetric 
systems and analytically some work has been done in the 
framework of the Thomas- Fermi approximation ^|,[7],|l3| . 
The negative scattering length case is mostly unexplored 
except some numerical results pj . Other approach to the 
dynamics of the condensate is the time dependent vari- 
ational technique Jl4| ] which assumes a fixed profile and 
computes the evolution of some parameters such as the 
width, etc., by variational techniques 

An important fact related to negative scattering length 
condensates is that stable solutions to Eq. (Q) exist only 
under certain conditions on the number of particles and 
the size of the trap |l4|-|r7|] . When those conditions are 
not satisfied the condensate is unstable and destroyed 
by the collapse phenomenon because the density \ip\ 2 in- 
creases up to a point where nonlinear losses (not included 
in Eq. (Q)) become dominant. So, to have a large stable 
negative scattering length condensate collapse must be 
avoided. Having larger condensates is important to get 
better experimental observations of BEC. The reason the 
critical number of particles that can be put in the con- 
densate without collapse is very small for current exper- 
imental parameters and thus, it is difficult to perform 
accurate measurements and to obtain experimental data 
of the condensation process. Other reason for the interest 
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of large condensates is their future practical applications 
(atom interferometers, atom clocks, etc.), where coherent 
atom clouds as large as possible are necessary. 

In this paper we concentrate on the analysis of negative 
scattering length condensates in cigar-shaped traps. We 
present a new class of soliton-like stable solutions which 
can be of interest in the applications of these coherent 
atom aggregates. 



II. DERIVATION OF THE MODEL EQUATIONS 

From now on we will study the solutions of Eq. ([!]) in 
cylindrically symmetric parabolic traps (||). More explic- 
itly we will consider cigar-like condensates, i.e., the case 
in which the trapping potential in s is much weaker than 
the trapping potential in p; mathematically A 2 < 1. 

Let us make the change of variables: r = vt^a^p = 
r, aos — z and Q = —8iraN/a where clq = \jTijmv is 
the size of the ground state solution of the noninteracting 
GPE with a harmonic potential of frequency v (except 
for a y/2 factor). Let us also define a new wavefunction 
as u(p, s, t) = ip(r, z, t)yJa^/N, then equation (pi) reads 



du 
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with the normalization condition J \u\ 2 df = 1. 

The solution of this nonlinear partial differential equa- 
tion is a challenging problem and no explicit solutions are 
known. However due to the different interaction scales 
involved in our particular problem it is possible to find 
approximate (but very accurate) analytic forms for the 
ground state solutions of Eq. (^) . A detailed analysis us- 
ing multiscale expansions is done in the Appendix. Here 
we will derive the ground state solution by simple phys- 
ical arguments. 

We will first assume that it is possible factor the solu- 
tion of Eq. (§) as 



u{p,s,r) = 0(p)£(s,r). 



(4) 



Then, <\> satisfies the following equation 

-p~(p = v p <t), (5) 



Eq. (||) is a well-known eigenvalue problem, the two di- 
mensional isotropic harmonic oscillator. Its ground state 
solution is 



MP) = e- pV2 . 



(6) 



Multiplying Eq. (||) by <j>* and integrating to eliminate 
the p dependence we find 
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where the additional factor 1/2 in the nonlinear term 
comes from the quotient J °° \(f)Q\ A pdp/ \tfio\ 2 pdp = 
1/2. 

Finally, let us make the change ip(s,r) — ^(s,r)e ll/pT 
to obtain 



(8) 



This equation is a 1-dimensional NLSE, that in the 
X z = case can be integrated by the inverse scatter- 
ing technique When X z = 0, Eq. (g) has stationary 
normalized single-soliton solutions of type 



(9) 
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From here and using the Galilean invariance of the ld- 
NLSE it is possible to find travelling soliton solutions 
that propagate without distortion. 

The width of the cloud in the s direction is related to 
the nonlinear coefficient through the relation 



< s z > L 
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This fact is remarkable and means that condensates with 
a small number of particles (Q is proportional to N) 
would be very long, while condensates with more parti- 
cles would be shorter. If the number of particles is large 
enough the condensate is unstable and collapse occurs. 

To justify approximation (Ji|) let us note that in the 
transverse direction the trapping potential and nonlin- 
ear force tend to compress the wavepacket competing 
against the linear dispersion effect provided by the ki- 
netic energy term. On the other hand the trapping force 
in the s-direction has been removed so that along that 
axis there is only a competition between the nonlinear 
attraction and the dispersion. When the main force on 
the transverse direction is the one caused by the trap- 
ping potential the approximation will be justified. To 
check it let us compare both potentials -fftrap = \p 2 and 
-ffseif-int. = for the soliton solution. Their ratio is 

a function q(p, s) given by 
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Since sech(a;) < 1/2, evidently when 32tt 2 /Q 2 3> 1 then 
q(p, s) 3> 1 except for very small values of p. When this 
condition is satisfied, the parabolic potential dominates 
over the self-interaction, and then the only effect of the 
nonlinear term on the transverse shape is to provide small 
shape corrections near the center of the trap, which is 
the place where the parabolic potential is lower and the 
nonlinear term more relevant. 

To see whether the soliton solutions really exist we 
have computed numerically the ground state solution of 
(^) for different values of X z . In the non interacting limit 
(small Q) the solution is given by 
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u( Pl s) = X 1 J i n~ 3 ^exp(-p 2 /2- 
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Decreasing A 2 and increasing Q preserving also the con- 
dition q(p, s) > 1 we should obtain the soliton solutions 
@ . To compute the ground state solutions we have used 
the steepest descent method described in f|| to minimize 
the Hamiltonian 



X = / 



(14) 
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over a discrete lattice, where the solution is defined. 
In Fig. 1 we plot sections of the fundamental state 
u(p, s) for different values of the trapping potential in s 
(parametrized by the value of A 2 ) and Q = 5. As the A 2 
value decreases the solutions get wider, but when A 2 = 
the solutions do not widen indefinitely. The profile is 
then very close to the one defined by Eq. (0), while the 
transverse profile is Gaussian as Eq. (H) predicted. For 
Q = 10 it is seen in Fig. 2 that the A 2 = solutions are 
not so close to the profiles predicted by Eq. (||). This is 
because when Q is large the approximation involved in 
the derivation of Eq. (J9|) is not valid and the nonlinear 
energy term is comparable to the transverse harmonic 
trapping energy. However, it is striking that even for 
this large Q case the differences between the numerically 
calculated profile and the solitonic profile are less than 
10% . In this case it is seen that compact solutions exist 
when the trapping in z is absent. 

On the other hand, when a strong trapping potential 
is applied in the s-direction (stronger than the nonlinear 
self-interaction term) the numerical solution is close to 
the exact Gaussian ground state solution (|l^), and the 
effect of the nonlinear term is only an enhanced com- 
pression of the solution near the center of the trap. This 
phenomenon is seen in Fig. 1 (plot with A 2 = 0.4). 

The existence of atomic solitons has been put forward 
in [p5[ in the context of the motion of an atom beam 
in the field of a travelling wave laser and in a similar 
context in p6[. In those papers however the trap effect 
was not considered and the validity of the transition to 
1-D equations was not studied. 



III. CONTROL OF THE CONDENSATE MOTION 

Thus, we have numerically established the existence 
of localized solutions and obtained their analytical form 
when the trapping in s is eliminated (provided the num- 
ber of particles N oc Q is small enough to avoid col- 
lapse). Also there exist travelling solutions of this type 
that could propagate without distortion. Now it is inter- 
esting to study the response of the center of mass of the 
condensate to an external potential to because it could 
allow to control the motion of the condensate. 

Let us then study the evolution of the center of mass of 
a condensate gover nedbyEq. (§) in an arbitrary external 
potential V{x,y,z) (the following results arc valid for any 
potential not only for parabolic traps). Defining 



and computing its time derivatives using Eq. (|l]) we find 
dX/dt =< V >, where V is the usual momentum opera- 
tor, V — ~ihV and 
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which is the Ehrenfest theorem of Quantum Mechanics. 
Eq. ( |l5| ) means that this theorem is still valid for the 
GPE so that the center of mass of the wavepacket behaves 
like a classical particle. It is possible to check the validity 
of ( |l~5| ) for more general NLS equations (i.e., more gen eral 
nonlinear terms) , a fact which is not well-known [|27| . 

This result implies that one could manipulate a con- 
densate by using an external potential as is known for the 
ld-NLSE p8fl . Joining this result with the previous one, 
i.e., the existence of localized solutions, we find a way to 
control the motion of a negative scattering length con- 
densate: just relax the trap in one direction and apply 
an external force along that axis, the condensate will re- 
spond by preserving its shape and moving like a classical 
particle. Of course the external force should be smoothly 
varying since the localized solutions have been derived in 
the limit where no forces are present |^9| . 

It is not strictly true that a condensate would respond 
as a whole to the external force. It is well known JlJ] 
that any initial data evolving following (||) decomposes 
into solitons and "radiation" . So, it should be simple to 
find experimentally these solitonic objects by just adia- 
batically relaxing the trap and applying an external po- 
tential. Doing so radiation (which in these context moans 
some free atoms) would be generated and some solitons 
obtained. Only if the initial data is already a soliton 
there would not be a breaking of the initial data into a 
soliton train plus radiation. 

Concerning the motion of a soliton wavepacket in a 
highly asymmetric trap (with A 2 small but different from 
zero) it must be said that there are no completely sta- 
tionary solutions as shown by |l9f| . However there would 
be a quasi-stationary solution with gaussian tails in the 
s — > 00 and near-sech profile in the s = region as re- 
cently proposed in pp| . If the number of particles is large 
enough so that the soliton size is small compared with 
the scale of variation of the potential it will be expected 
a smooth motion of the "soliton" towards the "bound- 
ary" . However, if the soliton size and trap size are of the 
same order of magnitude there will be a competition of 
both scales which can result in cuasiperiodic motion or 
even into chaotic motion as discussed in . 

Another "tool" to control the motion of the condensate 
could be a laser field as has been put forward in [p5 22 
but in those cases the interaction between the trasverse 
laser field and the atoms should be carefully considered. 
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IV. APPLICATION TO LITHIUM 
CONDENSATES 

Let us analyze the relevance of our results for the 
Lithium Bose-Einstein condensates P,|l2fl- Following 
Ref. |Q we will take as parameters a — —14.5 Aand 
the usual trapping potentials for the cigar-trap that are 
about v — 150 Hz corresponding to ao ~ 3/im. 

To ensure the validity of Eq. (|J) it is necessary that 
W s ^S> 1 and then we find that N ~ 300. However in 
Fig. 2 it is seen how even in the case Q = 10 (N ~ 900) 
the differences between the soliton profile and the real 
ground state are small. 

Another interesting limit corresponds to collapse. In 
principle one would expect that the cigar-like trap would 
allow a larger number of particles to be put before col- 
lapse occurs. The physical reason is that keeping free 
the condensate in one spatial direction collapse would 
not occur along that axis but through compression of 
the orthogonal (transverse) directions, which are smaller 
and thus 'feel' stronger interactions. This means that the 
system would behave in a two-dimensional like manner 
and then the collapse conditions should be less severe 

To test this hypothesis we have performed simula- 
tions of the largest Q value allowed using the same code 
as for the computation of the ground state. The up- 
per limit found for the cigar- like trap is Q — 17 cor- 
responding to N ~ 1500. This number is somewhat 
lower than the gaussian bound given in fl4j[ which is 
Q = 19.5 corresponding to N ~ 1710. These num- 
bers compare favourably with the spherically symmetric 
results. In that case the limit found using the steep- 
est descent method is Q = 13.7 corresponding to about 
N ~ 1200, again lower than the gaussian bound Q = 16.7 
and then N = 1460. So the cigar-like trap allows to in- 
crease the maximum number of particles by 25%. This 
is a small but significant increase. We have cheked by 
numerical solution of the gaussian equations of Ref. 
and numerical simulations of ([}]) that the cigar trap is 
the optimal one; i.e. there is no another parabolic trap 
configuration with better collapse-avoiding properties. 

As stated in the introduction one would like to increase 
the number of particles in the condensate as much as pos- 
sible. This would allow the control of a large coherent 
pulse of atoms. To do so one could try to use a higher 
order soliton po| . However, for those solutions the shape 
performs complicated (but periodic) oscillations and de- 
velops high spatial and temporal gradients that probably 
would rule out the approximation and induce collapse if 
the order of the soliton is large enough. It is also possible 
to generate a soliton train where the solitons have differ- 
ent global phases so that the interaction between them 
is repulsive. This idea should work for some situations 
allowing many particles to be put in the ground state and 
will be elaborated in future work. Other possibility is to 
use a non-Gaussian fundamental mode for the transverse 



solution as in |4|j31j]. However the question of the sta- 
bility of those solutions under general three-dimensional 
perturbations is not trivial and is the subject of current 
research. Finally, other possibilities proposed in the liter- 
ature could be of use here such as using two condensates 
[23 1 or the control of the value of the scattering length 
\24\. 



V. CONCLUSION 

We have found compact solutions of (||) that exist due 
to nonlinear effects even when the z trapping potential is 
absent. Joining this result with (|T^) we conclude that it is 
possible to control the motion of the condensate, which 
could propagate without distortion by using smoothly 
varying external potentials. Thus, the atom cloud could 
be manipulated very easily, e.g., with an atom guide. It 
is interesting and curious that this cigar-like packet could 
be transported in that rigid way behaving like a quasi- 
particle. This behavior is specific of negative scattering 
length condensates and an advantage over the positive 
scattering length ones, which tend to fill all the available 
space due to the repulsive atom-atom interaction. Addi- 
tionally we have pointed out that relaxing the trapping 
potential in one direction in current traps would allow to 
increase the number of particles that can be put into a 
negative scattering length condensate. 

We hope that this study will stimulate the experimen- 
tal efforts in performing BEC with negative scattering 
length and think that the soliton solutions here studied 
will be of practical applicability in Bose-Einstein conden- 
sate "engineering" . 
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APPENDIX: DERIVATION OF THE 1-D NLSE 
BY MULTIPLE SCALE ANALYSIS 

Here we will give the details of a more formal derivation 
of Eqs. (||) and (||) from (||) using multiple scale analysis 
[B2| . Let us first choose the 



u = e 1 / 2 Mo(r, t', x, y, z') + e 3 u\ + 
K = e 4 X z 
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Q = eQo 



(Ale) 



This scaling satisfies the desiderable property that the 
L 2 norm of u is conserved and that the potential in z 
is weaker than the nonlinear interaction (and the later 
wea ker t han the trasverse potential) when e — ► 0. Insert- 
ing (|Al|) into Eq. (§) we find 
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We now separe Eq. (A.2) in the different orders in e 
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where 



+ Eq. ( |A3a|) implies that the 

trasverse profile of uo is given by the isotropic two di- 
mensional harmonic oscillator equation and then uq can 
be ch oosen as uo = <j>(x, y)£,{z' , r')e l/T . Substituting into 
Eq. (A3a), multiplying by 4>* and integrating over the 
trasverse coordinates x, y we obtain 



1 d 
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This means that the longitudinal profile obeys the non- 
linear Schrodinger equation. In the X z = case the so- 
lutions can be found analitically as discussed in Sec. II. 
Joining the longitudinal and the trasverse solution and 
changing back to the nonscaled variables we find that the 
ground state solution has the form 



u(p,s,r) 



47T V 8tt 



(A5) 



at least to the first orde r in e oc Q. The corrections are 
given by Eqs. (A3c-A3i). It is easy to see that the equa- 
tions have solutions u\ = u 2 = 0, so the solution is de- 
termined at least to order e 11 / 2 by uq. This is the reason 
why the ground state sol ution is close to the approximate 
profile given by Eq. (A5) even in the nonperturbative re- 
gion as discussed in Sec. II. 
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FIG. 1. Sections of the ground state solution of Eq. (^) 
with Q = 5 and different values of \ z . From the innermost 
to the outermost curves the \ z parameter is \ z — 0.4, 0.2, 0.0 
The dash-dot line corresponds to the theoretical prediction 
for X z — given by Eq. (^) while the dashed line corresponds 
to the Gaussian solution given by Eq. (Q for \ z = 0.4. (a) 
s-section for p = 0. (b) p section for s = 0. 



FIG. 2. Sections of the ground state solution of Eq. (g) 
with Q = 10 and X z = 0. The dash-dot line corresponds to 
the theoretical prediction for X z = given by Eq. (^|) (a) 
s-section for p = 0. (b) p section for s — 0. 
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